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Abstract 

We find an exact spin liquid state without time reversal symmetry on the 
kagome lattice with odd number of electrons per unit cell and explicit 
wave functions for all eigenstates. We also obtain that all spin-spin 
correlations are zero except trivial cases. We then show that there are 
anyonic excitations in our model. Finally, we demonstrate the existence 
of ground state degeneracy and gapless edge states indicating nontrivial 
topological order in our model. We also label all eiginstates on torus by 
local and global string operators. 



PACS number(s): 75.10.Jm, 75.50.Mm 

Introduction. -There has been enormous interest recently in topological quantum phases 
[U [2] beyond the Landau paradigm [3] , which is based on the idea of symmetry breaking 
and local order parameters. The newly discovered phases have many exotic phenomena 
such as electron fractionalization [4 J and emergent gauge fields [5]. The traditional order 
parameter descriptions are either inapplicable or insufficient for these new phases. 

Spin liquid states in dimension two (2D) are examples of such exotic phases. In the spin 
liquid states, strong quantum fluctuation forbids existence of nonzero spin order parameter. 
Although many numerical and analytical evidences indicate the existence of such exotic 
state in 2D, exact solution is absent until very recently [HI EJ [101 [111 02] . These exact 
solutions provide us with a lot of intuition of the internal structure of spin liquid states. 

Spin liquid states that spontaneously break time reversal symmetry (TRS) and spacial 
inversion symmetry are called chiral spin states [13]. An exact chiral spin state is proposed 

* Corresponding author. Electronic address: slwan@ustc.edu.cn 



1 



(a) (b) 



Figure 1: (a) The kagome lattice is composed of two types of elementary plaquettes: tri- 
angles and hexagons. The dashed lines show the unit cell, (b) Elementary plaquettes of 
Kagome lattice. The letters x,y,z indicate the factors of F p . 

in Ref.[9] based on short-range anisotropic interaction [7\, but this model have even number 
of electrons per unit cell, and thus not deserve to be called "Mott insulators" or "spin liquid" 
states in the traditional sense. We consider it an interesting question to find whether there 
exist other exact spin liquid states with arbitrarily given symmetry properties and odd 
number of electrons per unit cell. 

In this letter we generalize the toric code model [8, 6j to the kagome lattice (Fig.l) and 
propose an exact spin liquid state that breaks TRS. This state is not chiral spin state in 
conventional sense because TRS is already broken in the Hamiltonian. The kagome lattice 
has odd number of spins per unit cell, and thus the proposed state is a true spin liquid state. 
Our model treat three spin components a x ' v ' z democratically, which is more symmetrical 
than the a x ' y model [8, 6j. We propose explicit wave functions for all the energy eigenstates 
and interpret these eigenstates as string-net [1] coherent states. We also compute spin-spin 
correlation functions and find that all of them are zero except the trivial cases of on-site 
correlations. After a short discussion of anyonic excitation in our model, we proceed to 
discuss the ground state degeneracy and gapless edge states indicating topological order. 
We also find all the eigenstates of the model on torus using string operators. 

The model. -Our spin-^ Hamiltonian on the kagome lattice is given as: 

H = -J2V P F P (1) 

p 

with summation over all elementary triangular and hexagonal plaquettes. We choose F p = 
a i a 2 a 3 on the triangles and F p = afa^o-^afa^ag on the hexagons (see Fig.l), where af are 
the well-known Pauli matrices. It follows from the property of Pauli matrices that F p = 1. 
Another crucial property of F p is that [F p , F q ] = for any pairs p and q. So all F p can have 
eigenvalues (±1) simultaneously. Since af change sign under time reversal transformation, 
the Hamiltonian breaks TRS explicitly. The solubility of our model does not depend on 
signs of V p , so we will choose all V p = V > for simplicity. 



Firstly we consider a planar region with N sites, where N is very large. The number 
of elementary triangle and hexagon are shown to be 2N/3 and N/3, respectively. So there 
are N plaquettes, and F p can label 2 N states. The number of physical spin states is also 
2^. From this counting argument we know that eigenvalues of F p label all the states on an 
infinite plane and the model is exactly solvable. 

We can write down the following wave unctions for all eigenstates: 



h/>A> = n( 1+ vWo> (2) 

p 

where p exhaust all elementary triangular and hexagonal plaquettes. is an arbitrary 
spin state which we choose as the state satisfying af = +1 for all sites. The irrelevant 
normalization factors have been omitted. Different choices of A p = ±1 label different states. 
The F p expectation in \ip\) is: 



<V>a|Va) ~ Xp (3) 

so |V>a) is an eigenstate of H with eigenvalue E = — J2 p X p V p . The ground state is the 
state with all \ p = 1. The above simple wave functions have transparent explanation in the 
physical language of string-net condensation [T]. We will demonstrate this in our model. 
A caution in order is that this wave function can completely determine the state only on 
an infinite plane. On systems with boundary like disk or closed like torus, there are other 
degrees of freedom missing. We will return to this question later. 

We define two types of strings T\ and T2 [U [T] . T\ strings walk on triangles while T2 
strings walk on hexagon (see Fig. 2). We define the string operators as Pt = IIiet a i f° r each 
string t. The choice of af depends on both the site % and the string type. For T\ strings, 
a a = a y ,a z ,a x at the three vertexes 1,2,3 of each triangle, respectively. For T2 strings, 
a a = a z ,a y ,a x , a z ,a y , a x at the six vertexes 1,2,3,4,5,6 of the hexagon, respectively. An 
important property of this definition is that T\ and T2 string operators involve different 
Pauli matrices at any sites i, from which we can readily obtain the simple commutative 
relations among Pt. For arbitrary T\ strings s,t and T2 strings u, v, we have [P s , Pf] = 0, 
= 0; [P S ,P U ] = if s and u intersect each other by even times, while {P S ,P U } = 
if s and u intersect each other by odd times. 

These strings are divided into two classes by their topology. A string is called open 
string if it has endpoints, otherwise is called closed string [lj. From the property of Pauli 
matrices we can find [Pt, F p ] = for any p for a closed string t, and it follows from this that 
[Pt,H] = 0. For open string s with two endpoints on plaquettes p, q, we have {P s ,F a } = 
if a = p or q, and [P s , F a ] = otherwise. So P s changes the sign of F p at its two endpoints. 
If we regard negative F p as vortex at p, then P s applied to the ground state will creates two 
vortexes at its endpoints and increase the energy by 2V. 

It is interesting that F p are also closed strings operators defined above. Actually, for a 
T\ string t\ surrounding a single hexagon p (such as string c in fig. 2), we have P^ = F p . 
Similar relations also hold for any T2 string surrounding a single triangle. So we conclude 
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Figure 2: Strings on the kagome lattice, a, b are open strings, while c, d are closed strings, 
a, c are T\ strings, while b, d are T2 strings. The letters x,y,z denote the specific Pauli 
matrices appearing in the string operators. 

that F p are special cases of closed string operators. This observation gives physical meaning 
to the above wave function. We expand the wave function: 

I^A) = (1 + E X p F p + E A P1 A P2 F P1 F P2 + • • -)|Vo) (4) 

P Pi ,P2 

each term is a closed string state. So \ip\) are coherent states of closed strings. 

We mention that there are also two types of strings in the models of Kitaev and Wen 
[U [8] , but strings in our model have the additional feature that T\ and T2 strings walk on 
triangles and hexagons, respectively. So their properties are different. One of the differences 
is that there are twice as many triangles as hexagons and thus more room for open T\ string 
endpoints than for T2 string endpoints. 

Exact spin correlation. -The spin-spin correlation functions can be extracted from the 
exact wave functions. Since wave functions of all eigenstats are known, we can obtained all 
the spin-spin correlations by direct computation: 

= " IK 1 + At^Ko^o) (5) 

s t 

Where s exhausts elementary plaquettes with the property {P s ,afo-j} = 0, and t exhausts 
elementary plaquettes with the property [P s ,afaj] = 0. We can check that such s exist 
except in the trivial cases when i = j and a = b. Since (1 + X S F S )(1 — X S F S ) = 0, we obtain 
the exact result: 

(V>aK4I^} = o (6) 

except the trivial cases (Vojof IV'a) = (V'aIIIV'a) = 1- This computation is analogous to 
the computation of spin-spin correlations in the Kitaev model |14j . where they are found 
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to be exactly zero beyond the nearest neighbor. The spin-spin correlations in our model 
are even more exotic. They are identically zero for two different sites including the case of 
nearest neighbors. We also mention that the spin correlations in the earlier models [TJ [B] 
can also be computed using explicit wavefucntions. Our method can also be generalized to 
computation of multi-spin correlations. Another feature is that the result is valid for all 
energy eigenstates, besides the ground state. This remarkable feature is also shared by the 
Kitaev's model on honeycomb lattice [71 114]. 

The zero spin correlations indicate very strong quantum fluctuation in our model. So 
the conventional idea of using spins as local order parameter will not give meaningful result 
here. But there is topological order pQ in this model manifesting in ground state degeneracy 
on closed manifold and gapless edge states at the boundary of open systems. 

Anyonic excitation. -The anyonic excitation in our model has similar properties to 
earlier models [U [7J [8] . The energy of an open string is not changed if we move the two 
endpoints far apart, i.e. the endpoints are deconfined. We call these endpoints T\ and T2 
vortexes. We also mention that these vortexes can be interpreted as Z2 vortex excitations 
of an emergent Z2 gauge field [U US] . 

The vortex hopping operators is very simple, they are the same matrices appearing 
in the definition of Pt- The statistical angles of T\ and T2 vortexes can be shown to be 
both zero by the statistical algebra method [16J. So the same type of vortexes see each 
other as bosons. The mutual statistical angle of T\ and T2 vortexes can be computed in a 
different way. Following Kitaev's work [6], we consider the physical process of moving T2 
vortex a around T\ vortex b along the path of string u (shown in Fig. 3). The initial state is 
= P s Pt\ip), where is a state with all vortexes far away from a, b. The final state is: 

|V>/> = P u\^) = PuPsPM) = -P 8 P t PuH) = (7) 

the result does not depends on detail shape of u. The only important thing is that b is 
enclosed by u. So the phase factor —1 is purely topological. The mutual statistical angle 
9 is defined by e 2td = — 1. So we have 6 = tt/2 which is half of fermion's. This anyonic 
mutual statistic is a signal of topological order in our model. 

Ground state degeneracy and gapless edge states. -The ground state degeneracy on a 
closed manifold often indicate topological order hidden in a phase [TTJ [fl [18]. To study 
the topological order of our model, we consider a periodic system, i.e. we put the system 
on a torus. For a m x n system (as Fig. 4), there are 2mn triangles, mn hexagons, and 
3mn sites (the sites at opposite boundaries are actually the same site and should not be 
re-counted). Now the mutual independency of F p is destroyed by the periodic boundary 
condition. Actually, one can readily check that JJp^triangie Fp = 1 an d Ilpg hexagon Fp = 1> 
so the number of independent F p are 3mn — 2. 

It has been mentioned that closed string operators Pt commute with H and thus do 
not change the energy eigenvalue. But does Pt have chance to transform a given energy 
eigenstate to another eigenstate that degenerate with it? For the cases when t is homotopi- 
cally trivial, i.e. t is able to be deformed to a point, the answer is no. Actually, when t 
is a homotopically trivial T\ string, it can be checked that Pt = Y\ p F p , with p exhausting 
the hexagons enclosed by t. So Pt cannot transform a given energy eigenstate to another. 
Similar result holds for T2 strings. 



5 



Figure 3: Statistic of open string endpoints. s is an open T2 string with one endpoint a and 
the other endpoint far away, while t is an open T\ string with one endpoint b and the other 
endpoint also far away, u is a closed string enclosing b. 

But there are also homotopically nontrivial strings called global strings on the torus. 
Consider four global strings t\, t2, £3 and £4 shown in Fig. 4. We have {P^jP^} = 
and {Pt 2 ,P^ 4 } = by their intersecting times. All other pairs of string operators are 
commutative. Since all these operators commute with F p , we can choose the common 
eigenstates of P^ , Pt 2 and F p as the basis of the Hilbert space. To find whether eigenvalues 
of Pt x , Pt 2 and F p can label all states completely, we just need to do a simple counting: 
All of Pt t , Pt 2 and F p have two eigenvalues ±1, so they can label 2 3mn ~ 2 x 2 x 2 = 2 3mn 
states. This is just the dimension of physical Hilbert space. So Pt ± , Pt 2 and F p are complete 
labels and we have found all energy eigenstates on the torus. From the above discussion 
we conclude that states with given F p has additional 4-fold degeneracy coming from two 
global string operators P%^ and Pt 2 . It follows from this that states with all F p = 1 are 
4-fold degenerate, i.e. the ground state degeneracy is 4. 

Furthermore, we can figure out how these eigenstates transform under Pt 3 and Pt 4 . Since 
{P t3 ,P tl } = and [Pt s ,Pt 2 ] = 0, we conclude that Pt 3 change the sign of Pt x but not of 
Pt 2 . By similar argument, we can also show that p 4 change the sign of Pt 2 but not of Pt ± . 
This argument also shows that both ±1 eigenvalues of Pt l and Pt 2 are accessible by some 
physical states. 

We mention that t\ and £2 can be chosen rather arbitrarily, with only requirements that 
they circle the torus in m and n directions respectively, and they are of the same string 
type. The physics here will not change with our choices because it is topological. 

We also find the existence of gapless edge states in our model, which is another man- 
ifestation of topological order. If we cut the torus open along the bold lines in Fig. 4, the 
dimension of Hilbert space is larger by a factor 2 L than the case of torus, where L ~ m + n 
when m, n are large. But the number of F p labeling energis is unchanged when we cut the 
torus open. So there is a large degeneracy proportional to 2 L , which we can interpret as 
various edge excitations with zero energy [HE]. 

Conclusion. -An exact spin liquid state has been given on the kagome lattice based on 
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Figure 4: Kagome lattice on torus, m = n = 3 in this figure. The effective region is within 
the dark bold lines. The four copies of the point ll A" should be regarded as one and the 
same point because of periodic boundary condition. t%, t 2 , t% and t± are all closed strings 
on the torus. t\, t 2 are T\ strings, while £3, £4 are T 2 strings. None of them can be deformed 
to a point and thus they are called global strings. 

multi-spin interactions. All spin-spin correlations are found to be zero except trivial cases. 
There are also anyonic excitations and ground state degeneracy indicating nontrivial topo- 
logical order. We think that this exact spin liquid state with odd number of spins per unit 
cell might be helpful to improve our understanding of 2D spin liquid states. Furthermore, 
the existence of such exact states suggests that there maybe exist other interesting states 
based on multi-spin interactions. 

The work was supported by National Natural Science Foundation of China under Grant 
No. 10675108. 
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